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Abstract 

A global view is given upon the study of collapsing shear-free per- 
fect fluid spheres with heat flow. We apply a compact formalism, 
which simplifies the isotropy condition and the condition for confor- 
mal flatness. The formulas for the characteristics of the model are 
straight and tractable. This formalism also presents the simplest pos- 
sible version of the main junction condition, demonstrated explicitly 
for conformally flat and geodesic solutions. It gives the right functions 
to disentangle this condition into well known differential equations like 
those of Abel, Riccati, Bernoulli and the linear one. It yields an alter- 
native derivation of the general solution with functionally dependent 
metric components. We bring together the results for static and time- 
dependent models to describe six generating functions of the general 
solution to the isotropy equation. Their common features and rela- 
tions between them are elucidated. A general formula for separable 
solutions is given, incorporating collapse to a black hole or to a naked 
singularity. 

1 Introduction 

Spherically symmetric radiating spacetimes are important in astrophysics for 
modelling radiating stars and in cosmology. Gravitational collapse is a highly 
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dissipative process, required to account for the enormous binding energy of 
the resulting object Pp. In the diffusion approximation this is described by 
the heat flux. It allows to join the interior solution to the Vaidya shining 
star exterior [2]. Radiating models are necessary in cosmology to describe 
the formation of structure, evolution of voids and the study of singular isties 

i- 

Shear-free perfect fluids with heat flux are often studied in order to sim- 
plify the calculations and allow realistic analytic solutions. Two of their 
advantages are that there are just two metric components and their space 
evolution is governed by the isotropy condition, which is an ordinary sec- 
ond order linear differential equation in the radial variable [3], [I]. In the 
presence of heat flux the only non-trivial non-diagonal component of the 
Einstein equations becomes an expression for it. The vanishing of the heat 
flux implies a severe constraint, which transforms the isotropy condition into 
a non-linear and highly complicated differential equation with few explicit 
solutions [3] , [5] , [6] , [7] . This situation persists even when anisotropy of the 
pressure is allowed both for shear-free and geodesic fluids [8], [9]. 

Many analytic solutions of the isotropy condition have been found. It 
has been written in a compact form already in 1948 [5]. Published in an 
obscure journal, this paper has been discussed nevertheless in some mono- 
graphs [3], [TJ. In spite of this, researchers prefer to work directly with the 
metric coefficients, which complicates the investigation. The isotropy condi- 
tion is the same for static perfect fluid models and for time-dependent ones. 
In the static case both the heat flux and the off-diagonal component of the 
Einstein tensor vanish identically, so there is no additional constraint on the 
isotropy condition. Time dependence is obtained by promoting the integra- 
tion constants into arbitrary functions of time. As a result, there are two 
groups of authors - the static (S) and the dynamical (D) group. Strangely 
enough, almost no interaction exists between them. This problem becomes 
especially annoying when the general solution of the isotropy equation is dis- 
cussed. Recently, such a solution was proposed by D authors [10], using the 
Lie symmetries method. It was treated as a class of solutions at first, but 
later this statement was corrected [TTJ. However, generating functions have 
been found by the S group as early as 1971 [12], [13]. Interestingly, the S 
authors were not aware of previous research inside their group, so more such 
functions appeared during the years [Il],[l5],[T6],[T7J,[T8] (in the last refer- 
ence a connection is made between it and the previous one). Of course, only 
non-vanishing in the static limit characteristics of the models were studied, 
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like energy density, pressure and the mass function. 

While in the static case the main junction condition to the Schwarzschild 
solution is the vanishing of the pressure on the surface, in the dynamical case 
the correct joining was found in 1985 [2]. Like the isotropy condition, this 
is an ordinary differential equation. The difference is that it has only time 
derivatives and is non-linear. Once again, mainly its formulation in terms 
of the metric coefficients was considered [IS] , [20] , [21] • During the solving 
process it has been found that some combinations simplify the computations. 

One of the important questions of gravitational collapse is whether it ends 
in a black hole or a naked singularity. A model with separable metric was 
thoroughly studied [22] , [23] , [21] and it was found that a black hole forms 
at the end. However, there exists a simple solution of the same junction 
condition, when horizon never appears and the fate of the collapsing matter 
is a naked singularity [25J. Further, conformally flat and geodesic models 
were studied extensively, but the fact that the latter are a subclass of the 
first is not universally known. 

In the present paper we address all these questions. In Sect. 2 the field 
equations are given as well as the definitions of important characteristics of 
the fluid spheres like energy density, pressure, heat flux, expansion scalar, 
mass function, the second Weyl invariant and the distribution of the causal 
temperature. In Sect. 3 the junction conditions between the interior and the 
exterior are presented and the definitions of quantities that lie or depend on 
the surface are written. These are the surface luminosity, redshift and tem- 
perature, the luminosity at infinity and the total energy radiated during the 
collapse. It is stressed that the energy is stored in the integration functions. 
The condition for the formation of a horizon and consequently a black hole 
is also mentioned. Sect. 4 contains the compact formulation of the isotropy 
equation (LG formalism) and its general solution, obtained in six different 
ways. The relations between the different generating functions are elucidated 
and the characteristics of the model are expressed through some of them. In 
the next few sections several classes of solutions are derived utilizing the LG 
formalism. These are separable solutions in Sect. 5, conformally flat and 
geodesic solutions in Sect. 6 and solutions with functional dependence be- 
tween the metric components in Sect. 7. For the sake of completeness we 
present chronologically the remaining solutions known to us in Sect. 8. The 
last section contains some conclusions. 



3 



2 Field equations 

The collapse of a shear-free perfect fluid sphere is described by the following 
metric in isotropic comoving coordinates 

ds 2 = -A 2 dt 2 + B 2 (dr 2 + rW) , d^l 2 = d6 2 + sin 2 Odtp 2 , (1) 

where A and B depend on r and t. The energy-momentum tensor for a fluid 
undergoing dissipation in the form of heat flow reads [4J,[22J 

T ik = {^ + p) UiU k + pg ik + qiu k + q k u h (2) 

where /i is the energy density of the fluid, p is the isotropic pressure, Wj is 
the four-velocity and is the radial heat flux vector, which is orthogonal to 
Ui. In comoving coordinates 

u^A-Hl g* = (0,g,0,0). (3) 

The non-trivial Einstein equations are 

3B 2 1 (2B" B' 2 AB'\ 

1 / 2B B 2 2AB\ 1 (B 12 2A'B' 2 A' 2B'\ . , 
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^=T*(tb)- ® 

Here the dot and the prime stand for time and radial derivatives respectively 
The rate of collapse G is given by [19] 



Let us use the variable u = r 2 and organize the time derivative terms in 
terms of 0. The Einstein equations become 

torn Q2 4 ( Wu uB l, 2uB uu\ , Q . 



8np = — — 



e 2 

3 



26 4 
3A + B 2 



uB„ 



B 



— + 
B 



2A 
A 



IBl 



+ 



8vrg : 
2A U B U 



4r9 M 
3B 2 

A 



B,, 



Aqt 



0. 



~B 



(10) 

(11) 
(12) 



B 2 ' AB A B 
The last equation is the difference between Eq (5) and Eq (6) and represents 
the isotropy of pressure. It contains no time derivatives and is an ordinary 
second-order differential equation for A, B. Thus we have the freedom to 
choose arbitrarily one of them and solve for the other. Time dependence 
arises when the integration constants Cf are promoted to integration func- 
tions Ci(t). Eqs (8-11) become expressions for Q,p,,p and q. The effective 
adiabatic index of the fluid T =dlnp/dln/z also can be found. 

When the heat flow q vanishes Eq (11) becomes another condition on 
A, B and the combination with the isotropy condition leads to a non-linear 
equation with few solutions [3], [5], [6], [7]. This situation remains also in the 
anisotropic case for shear- free or geodesic collapsing spheres [8], [9]. 

The static case follows when the functions Cj (t) become constants again 
and = = A. Then Eq (11) yields simply q = but the isotropy condition 
remains the same. Thus to every dynamical model corresponds a static one. 

An important characteristic in the general case is the mass function m of 
the fluid ball [19] 
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The conformal tensor has one essential component, which is given in an 
invariant way by the second Weyl invariant ^2- It can be determined from 
the following formula, holding for anisotropic fluid spheres with heat flow [8] 

An 

(H+Pt-Pr)-*2- (14) 



m 

In the isotropic and shear-free case R = rB and we have 

87r/i m 



Eqs (9,13) yield 
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The terms with time derivatives cancel and the equation for conformally flat 
solutions ^2 — is an ordinary differential equation in u, like the isotropy 
condition (12). Formula (14) has been derived originally in terms of the 
electric part of the Weyl tensor E [26J . Comparing both formulas one comes 
to the conclusion that E = 3^2- 

Another characteristic of the collapsing sphere is the temperature distri- 
bution T among its volume. Unlike the previous quantities, it is determined 
by a non-linear differential causal transport equation. For the present metric 
it becomes [20] 

T(qBy + qAB = -^l, (17) 

where k is the thermal conductivity and r is the relaxational time-scale which 
gives rise to the causal behaviour of the theory. Both of them depend on the 
temperature in general. A physically reasonable choice is the transportation 
of energy by massless particles. Then Eq (17) becomes 

jX' + qA 4 - a B 2 X^ 4 + P(qB)- A 3 B = 0, X = (AT)\ (18) 

where a, (3, a are constants. In the non-causal case — it is easily solved. 
Integrable causal cases are a = 0;2;4 [1] , [20] , [2T] , [28J , [29J . Explicit solutions 
do not alter the structure of this equation. 



3 Junction conditions 

The collapsing fluid lies within the sphere £ defined by r = r^. The fluid 
is radiating, hence, the exterior is not vacuum, but the outgoing Vaidya 
spacetime with the metric [I] 

ds 2 = - (l - ™M\ dv * _ 2dvdp + p 2 dn 2 . (19) 

The junction conditions represent the continuity of the first and second fun- 
damental forms on S. This results in the relations [2] 

p = qB, M (v) = m(r,t) , p{v) = rB (20) 

which hold on the surface S. Here M is the total mass of the radiating 
sphere, ps is its radius as seen from outside and v is the time of the distant 
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observer. The first relation is a differential equation containing derivatives 
with respect to t only. Replacing Eqs (10,11) and setting r = r s we get 



^ 26 4r6„ . , 

T + — = Psi (21) 

where p s is the pressure of the corresponding static model. 

Some of the characteristics of the model are defined on its surface. Such 
are the surface luminosity A s and the redshift z% [4|,[30j 

3/2, /. 2uB„ rB<d\~ 1 



A s = -uT (BQ U ) S , ^={1 + -^ + — -1. (22) 

The exterior time v is related to the interior one t as follows 

v = f (1 + z s ) A E dt. (23) 

The surface temperature of the star is 

T ^ = i^| = (e^) E ' (24) 

where 5 is some constant. The total luminosity for an observer at rest at 
infinity reads 

_dM 2 f / 2uBj, rBey 

cfc (l + z E ) 3 V 5 3 / E ' 1 ; 

The total energy radiated during the collapse of the fluid sphere follows 
from the previous equation 



(•Ve 

= / A^dv = M (v b ) - M (v e ) , (26) 



where v\> (v e ) is the exterior time of the collapse's start (end). These cor- 
respond to t b (t e ) according to Eq (23). Eqs (13,20) show that the radi- 
ated energy results from the change in the integration functions ACj = 
Ci (h) — Ci (t e ). Thus the energy that a static model can give out is stored in 
its constants, which are animated during the collapse and evolve under the 
single condition given by Eq (21). 
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A separable solution was discussed in detail [22] , [23] , [23] . There tb = — oo 
and the exterior solution at that moment is the static exterior Schwarzschild 
solution in isotropic coordinates 




Here M is the constant total mass. The sphere collapses until a black hole 
is formed at t e = tsH- This happens when the coefficient goo m the exterior 
metric (19) vanishes at the surface £ and a horizon appears. Eq (20) shows 
that the relation 2m = rB should be satisfied there. It may be written with 
the help of Eq (13) as follows 

The first multiplier vanishes in order to satisfy this relation. Then Eqs (22, 
25) lead to the blowing up of the redshift and vanishing of the luminosity at 
infinity. 

Another separable solution with t& = —00, however, never develops a 
horizon [25J and the collapse proceeds till all the mass is burnt out, namely 
M (te) = 0, which maximizes E^. The final state can also be any static 
model for anisotropic fluids with vanishing radial pressure p r = and no 
heat flow [5T] . In this case Eq (21) is trivially satisfied. 



4 Six generating functions 

All of the previous formulas are based on the solution of Eq (12). The 
search for its solutions spans an interval of 63 years, starting from 1948 
when Kustaanheimo and Qvist wrote it in a very compact form [3], [5], [7]. 
Introducing instead of A and B the potentials L — 1/B and G = A/B it 
becomes 

2GL UU = LG UU , B = l/L, A — G/L. (29) 

This is a linear second-order differential equation. One can choose an ansatz 
for G and solve for L or vice versa. However, a general solution is hard to 
find. Choosing the function K = L u /L transforms Eq (29) into a Riccati 
equation for K 

K u + K 2 -%^ = 0, (30) 
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which is first order, but still a general solution for any G is not known. 

This difficulty may be overcome if we choose one of the potentials as a 
function of AB or its w-derivative. For example let us take A and W = 1/AB. 
Then 

L = AW, G = A 2 W (31) 

and Eq (29) becomes 

T^yiF' B = l ' AW - ( 32 ) 
It is readily integrable when an arbitrary W is given [TO] . 



A = A (t)ex V ± J sji^du. (33) 

The function of integration Aq (t) may be removed from A by a time change 
but it remains in the expression for B. In the above reference the solution 
was presented as a result of Lie symmetries analyses, together with other 
classes of solutions. Later it was emphasized that this is a general solution 
[TT] and in this way the potential of Msomi, Govinder and Maharaj W is a 
generating function for the metric A, B and all of the above characteristics 
of the model can be expressed through it. In addition, W together with Eqs 
(32-33) comprises the general solution of the isotropy condition Eq (29). The 
potential W = W' 1 = AB was also discussed [JO] with similar expressions 
for A, B. 

The story does not begin here, however. As emphasized, the isotropy 
condition holds also for the static case, which was studied extensively in the 
past. Unfortunately, most of the authors worked in the so-called curvature 
coordinates. Yet there is some amount of papers in isotropic coordinates. 
Quite a few concrete solutions have been found and beside them five other 
generating functions. As we pointed out in the introduction, the work of this 
'static' group of authors is completely unknown to the 'dynamical' group, 
which studied the time-dependent metric. The opposite is also true. One of 
the purposes of the present paper is to present both the static and dynamical 
results together, so that the future efforts may be united and rediscoveries 
avoided. We shall derive the generating functions from one another. Amaz- 
ingly, their authors were unaware of the work of each other with one small 
exception. 
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Thus, let us take instead of W the potential 



U = -^ = 2(lnAB) u , (34) 
U U = \U 2 -H\ # 2 = ^f- (35) 



In this form the isotropy condition was given by Kuchowicz [12] , [13] , [32] . For 
U this is a Riccati equation, but for H is an algebraic one and for A is a 
simple linear equation. The metric is given by 

A = exp^J Hdu, B = B (t) exp -J (U - H) du, (36) 

where Bq is a function of integration. Thus U is another generating function. 
The change H 2 = UJ turns the Riccati equation for U into a Bernoulli 
equation, which is integrable, or into a simple algebraic equation for J 

U u = l -U 2 - JU. (37) 

Here one can choose either U or J as a generating function. 

Let us replace next the potentials U, H by /, g according to the following 
expressions 

U = T H = T (38) 

so that 

7 = 2 ^- » = wxk u (39) 

Then the isotropy condition becomes simply 

fu = g 2 - \ (40) 

and one can choose either / or g as a generating function. The other one 
follows immediately, while the metric is given by 

A = exp - J jdu, B = fio exp i J 9 du. (41) 

This third generating function was found by Goldman [14], who also gave 
some particular solutions. His work was corrected and further developed by 
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Knutsen [TH]. The latter author expressed the characteristics of the static 
model in term of the potentials, studied the energy conditions and proposed 
another particular solution. 

Next, let us take the potential 



$ = lnA5=/0. (42) 



Then Eq (40) yields 



where e — ±1 and the metric is expressed through $ 



9 = , (43) 



A = expe— J - <5> uu du, B = B exp I $ - e— j 

(44) 

These are essentially the expressions of Ref . [18] and $ is the fourth generating 
potential, proposed by Lake. 

The fifth one Z was introduced by Rahman and Visser [17] . It is obtained 
from the relation 

$ = 2 / ^ rfw (45) 



1-uZ 

and is connected to the Goldman-Knutsen potentials by the equations 

2f = --u, g 2 = -— . (46) 
Inserting them in Eq (41) we get for the metric 

1 f V — Z u _,_ n „ ,_1 . f Z 



A = exp ±—= / -du, B = BqA' 1 exp / du. (47) 

a/2 J 1 — uZ J 1 — 



Finally, we should mention the potentials introduced by Stewart [15] 

P = 2r (InAB) 5 = 2r f In 4") , (48) 

V B J u 

which satisfy the equation 

2rP u -j-^P 2 + SP+ ^S 2 = 0. (49) 
11 



It can be shown that 



-, S=(2g-l)P 



(50) 



and then Eq (49) transforms into Eq (40). In this way P is the sixth (and 
the last known to us) generating function. 

Eq (40) shows that to every function / correspond two functions ±g, 
Suppose we have a solution A ly Bi with potentials f\,g\. Then there is 
another solution with f 2 = f\ and g 2 = — g\- Eq (41) yields 

A 2 = A{\ B 2 = B l A\. (51) 

This is exactly the Buchdahl theorem |33j in the spherically symmetric case. 
Now let us give the characteristics of the fluid model in terms of L and 

G: 



ds 2 
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-G 2 dt 2 + dr 2 + rW 
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(56) 
(57) 
(58) 
(59) 



It is clear that G appears only through O, except in the pressure. Note also 
the simple formula for the second Weyl invariant. 

Finally, let us express the characteristics in terms of a generating function. 
The most convenient are the Goldman-Knutsen potentials. A and B are 
found from Eq (41). The rest are 



2A 



du, 



(60) 
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(61) 

(62) 

(63) 
(64) 

(65) 



while g, As and are given by Eqs (11,22,24) respectively. In the junction 
condition (21) one should put 



f 2 B 2 



2f + u(l-g 



(66) 



computed at the surface S. Here g is given by Eq (41) and / (r,t) is an 
arbitrary function of u and a number of functions Cj (t). When the latter 
and B (t) are constant, q, 0, As, Aoo, Ts, T, vanish, so that the model 
becomes static. Then the formulas for the energy density, the pressure and 
the mass coincide with those of Knutsen |16|. 



5 Separable solutions 

Separable solutions have the following metric 

A = j (t) a (u) , B = h (t) (3 (u) . (67) 

The function j (t) may be set to 1 by a time change. Then Eq (12) shows 
that a and f3 should satisfy the isotropy condition as a static metric. Thus 
the generating functions are static, while h (t) = B (t). Replacing the above 
metric in the expressions for the various characteristics we obtain 

3h 4ra u h 
9 = ^' 8nq = -^h^ (68) 
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A — PH. r --^fel (69) 

3/i 2 87r/i s A 2 + 2M 87rp 5 

87F/i= a^ + ^' 87rp= -^^ + l^' (70) 

m = ^Lhh 2 + hm s , * a =^, (71) 
2cr /r 

2S U/3 + 2Ma/3 u -r/3 2 /iy s ' 

where quantities with an index s correspond to the static model with metric 

a, (3. 

We suppose that the boundary of the static model is also at r s , so there 
p s = and Eq (21) becomes 

2hh + h 2 -2ah = } a = 2 {^j^j ■ (74) 

This ordinary second-order differential equation governs the behaviour of h. 
Its first integral reads 

h=A= (aVh - b) , (75) 
where b is an integration constant. Integrating once more gives 

h b i— b 2 , i— b , , . 

t-t = — + —Vh+— In Vh--. 76 

2a a z a 6 a 

Here to is a second integration constant. When b ^ one can enforce the 
equality b = a by absorbing a constant in /3 [22],[23J,[24J,[28],[22]. This 
solution was analysed in great detail and leads to the formation of a black 
hole. 

When 6 = the solution is h — 2a(t — t ) [25J. In this reference this 
particular solution was given as a simple ansatz satisfying Eq (74), the general 
solution being unavailable, according to the authors. Here we have derived it 



14 



from the general formalism. It is easy to see that m and rB are proportional 
to t — to, hence, their ratio is constant. At X 



/2m\ 

\7b) 



E 




(77) 



This certainly can be made less than unity by choosing the constants of 
the arbitrary static solution. Going back to the arguments that lead to 
Eq (28), one comes to the conclusion that no horizon is formed during the 
process of collapse. Its end is marked by t e = to when the mass burns out 
completely and vanishes. The energy accumulated during the collapse is 
radiated at the same rate. Both luminosities, the surface temperature and 
redshift are constant, while fi,p and \J/ 2 diverge as (t — t )~ 2 , q ~ (t — t )~ 3 
and ~ (t — to) . This may be an indication for the formation of a naked 
singularity. In Ref. [25J the special solution 

a = 1 + cu, (3 = 1 (78) 

was considered, which in addition is conformally flat. Solutions with no 
horizon appear also in higher dimensional spacetimes [35].|36j. 



6 Conformally flat solutions 

These solutions have ^2 = and a look at Eq (56) shows that the LG 
formalism is the most appropriate for their study. Eqs (29,56) yield 

Luu = 0, G uu = 0. (79) 

Integration produces four integration functions. In the general case they are 
independent and Eq (52) indicates that in G one of them may be set to unity. 
Hence 

L = Cx(t)u + C 2 (t), G = u + C 3 (t) . (80) 

The metric and the combination AB, which is the basis of the six generating 
functions, become 

u + C 3 _ 1 u + C 3 

A -c^c^ B -c^Tc- 2 > AB ~ ( Cl u + c 2 r (81) 
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The characteristics of the fluid model depend on the three functions Cj (t). 
It is much simpler to use C\, L and G instead. One should note that 



L U = C U G u = l, 6 u = -^ + ^. (82) 

We can use Eqs (52-59) for the characteristics of the model, inserting in them 
the above formulas. Eq (56) yields 

m = -£r- (83) 



which follows also from Eq (15) or the square of the conformal tensor [19J. 
The time evolution of Ci,L and G is governed by the junction condition Eq 
(21) on S, where p s is given by Eq (55) with = 0. Dropping for a while 
the index E, we obtain after a lengthy calculation 



2LLG + [3L Z - 2LL) G + ArLLG - ArdLC - 4 (L - 2nd) LG 2 + 

+ (8C1L - UuCf) G 3 = 0. (84) 
With respect to G this is an Abel equation of the first kind 

AiG + A 2 G + A 3 G 2 + A A G 3 = 0. (85) 

It is not soluble analytically in general, but in some degenerate cases one 
obtains [21 J an algebraic equation (when A\ = 0) or Bernoulli equations 
(when A 3 or A4 vanishes), which are integrable. 
With respect to C\ this is a Riccati equation 

RxCx + R 2 C 1 + R 3 Cf + R A = 0, (86) 

which may be transformed into a inhomogeneous second-order linear equa- 
tion. The general solution may be found if a particular one is known. 

Finally, in order to elucidate the character of Eq (84) with respect to L 
we make the replacement 

L = T 2 , (87) 

which gives 

Gl -(G + 2rG) 1 - G 2 l - rG 2 l 3 C 1 + 2G 2 (u + G) l 3 C x - 
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- 3uG s rCf = 0. (88) 

One can make this equation linear and homogenous in I in two ways. First, 
we simply put C\ = [19], [20]. When G is constant there are three classes 
of solutions. In these references G was taken in the form G = C%u + 1, 
which leads to more involved coefficients of the equations. They hold for u in 
general, but are needed for w s only. Another possibility is to make L ~ C\ 
This happens when C*2 = kC\ with k some coefficient. The case of constant 
G was studied pp. In fact, in these two cases A and B become separable and 
Eq (88) is the analogue of the integrable Eq (74). Some other conformally 
flat and separable solutions with a single, linear in time, integration function 
have been discussed [25], [57]. 

A well known example of a static conformally flat solution is the interior 
Schwarzschild metric. In isotropic coordinates it looks like [TTJ, [38J 

I _|_ 1+C2 " i 

A = B = — (89) 

1 + u/ci 1 + u/ci 

and has two constants. The limit c 2 — > oo leads to the Einstein universe, 
while c 2 = —1/2 yields the De Sitter universe [17] . The Z potentials of 
these solutions were given in the last reference. As pointed out in the begin- 
ning, one can 'animate' these classical static models by making the constants 
time-dependent. A model where they, in addition, depend on each other 
was given by Kramer [39] and studied later [3Q]. The junction condition 
becomes a complicated nonlinear second-order differential equation, which 
surprisingly may be solved in terms of a special function. Due to their sim- 
ple structure, conformally flat solutions were among the first to be discovered 
0,@U,02],03,H>P!>0S]- The subclass with G = 1 was studied too [37]. 

An important class of solutions are the geodesic or non-accelerating ones, 
which have A = 1. This leads to G = L and the metric becomes 

ds 2 = -dt 2 + L~ 2 (dr 2 + r 2 dn 2 ) . (90) 

Eq (29) gives L uu = and Eq (56) shows that the geodesic solutions are a 
subclass of the conformally flat solutions. We have for the expansion 
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The characteristics of the model are obtained when these formulas are in- 
serted in Eqs (53-59). For example, the pressure reads 

B 2 2<~) 

8np= ALd + AuCf. (92) 

3 3 

The condition A = 1 simplifies the generating functions of geodesic models. 
Thus W = L, H = 0, g = 0, $ = - InL, Z = const and P = -S. 
The junction condition Eq (21) becomes 

- 4rL 2 C l + AL (rL + L 2 ) C x - AuL 2 C 2 = 2LL - 5L 2 , (93) 

taken on the surface S. It differs substantially from Eq (84), although it still 
represents a Riccati equation for C\ jH]. When some of the coefficients are 
made to vanish, a Bernoulli or a soluble Riccati equation follows. The ansatz 

L = L (t + t )\ (94) 

where Lo,to-n are constants, transforms Eq (93) into the linear second-order 
equation of the confluent hypergeometric function. Elementary solutions 
were found for n = 0; —2/3; —2 [48]. This method also reproduces the solu- 
tion with 

C\ / 1 — C20sexps\ 2 /6t \ 1//3 



B = — — ZJ ^ s 2 , 8 = - , (95) 

2c 2 V 1 — uc 3 exp s) V ci / 



where c, are constants. It has been discussed extensively in the past 
[49] , [50] , [51] . Finally, it should be mentioned that the study of collapsing 
shear-free perfect fluid models with heat flow began with the well-known 
Robert son- Walker cosmological model by promoting its constant A; to a func- 
tionA:(t) 0,0U, [ISl&P]. 

7 Functional dependence between A and B 

A class of solutions to the isotropy condition has the functional dependence 
A (B). The latter is equivalent to G (L). Then Eq (29) becomes 

2GL UU = GllLL^ + GlLL uu , (96) 

which may be written as 

^Zg-L u = ^. (97) 



L 



'L 
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Integrating once we obtain 

L u = d (t) exp / J* LL dL, (98) 
J —-G L 

where G\ (t) is an integration function. A second integration gives 

J exp (- J dLj dL = G x (t) u + C 2 [t) , (99) 

Ci (t) being a second integration function. Choosing an explicit G (L) one 
obtains after the integrations an explicit or implicit expression for L(t,u). 
The result for the variables A, L is similar 

/ 6XP (~ / 2A A-LA L L L dL ) dL = Cl {t) U + ° 2 {t) ■ (1 ° 0) 

This formula was found long ago [13] and considered to be the general solu- 
tion of the isotropy equation. Formally this is true, because time plays the 
role of a parameter and not a variable in it. Therefore, effectively, A = A (u) 
and B = B (u). Inverting the second equality and replacing it in the first, one 
finds that for every solution indirectly A = A (L). Thus L may be considered 
as another generating function. However, the inversion and the double inte- 
gration in Eqs (99,100) make the procedure rather cumbersome and implicit. 
Formula (100) was later rediscovered by the Lie symmetry method [10J and 



a few examples were presented for illustration in both references. In the LG 
formalism one of them is given by G = L 3 . A brief calculation yields from 
Eq (99) 

L=[C 1 (t)u + C 2 (t)] 1/7 . (101) 

When A is a function of B (or L) such is the base for the generating functions 
AB. Then Eq (39) shows that the Goldman-Knutsen potential g = g(L). 
However, / depends on L and L u in the general case. 

A model of the same type was found in an attempt to generalize the 
exterior Schwarzschild solution (27) [51], [IS]- It has 

^4 = ^4 B = B (t) (1 + Ff , F= Cl (102) 
1 + ? (1 + c 2 u) 1 

Bayin [56] studied static solutions in isotropic coordinates of the type 

A = A (p- c \ B = B Q< p C2 (103) 
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with 

= c 3 e C4U , (j) 1 - = c 3 u + c 4 , (104) 

where c = c(ci,C2) in a specified way. They also have A and B directly 
dependent on each other. 

8 Other solutions 

As stated before, static and time- dependent solutions are on equal footing 
with respect to the isotropy condition. We present only one of the metric 
functions in most cases, preserving the original notation. The constants 
below are understood. Chronologically, the first solutions were given by 
Narlikar et al in 1943 [Z] , [13] , [5Z] • They are static and have 

Lj = C ir 1+n/2 + C 2 r l - n/ \ L n = C 3 r- k/2 , (105) 

where n > 0, > k > —2. There are three cases of A for each L, including the 
interior Schwarzschild metric (SIM). Nariai [I],[5B] found five static solutions, 
one of them coinciding with Lj from above. The rest are 

L 2 = (a + bu) a , L 2 n = u (a + 6 ln r) 2 , L 2 in = a cos (b + cu) , 

Aj v = {a + bu) a . (106) 

In 1968 Strobel [3].|41j gave a list of cases of L uu /L for which a solution 
of Eq (29) can be found in the handbooks on differential equations and 
two explicit solutions with L uu = 0. This reference remained unnoticed. 
Kuchowicz [Z],[l2],[T3],[32j,[59j gave a host of new static solutions, using his 
generating function. In view of their great number we direct the reader to 
the original references. Goldman [13] studied three explicit static examples 
of his potential g 

gi = — , gn = —= coth (a — bu) , gm = cosh (a + bu) . (107) 

1 — bu v2 

Stewart [7],[15] applied the Buchdahl theorem to SIM to obtain a new static 
solution with 

A 2 = c(l + au) 2 (l-buy 2 . (108) 
Sanyal and Ray [43] gave their Case 1 dynamical solution 

A = C(t)u + D(t) (109) 
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as a complementary one to their general solution (100). Modak [H] proposed 
a time- dependent metric, which coincides with the animated fourth Nariai 
solution. Pant and Sah [60J studied in detail the static model with 

A = A B = — — , 6{u) = - 1/2 , 110 

1 + M 1 + u/a 2 (l + bu/a 2 ) 1 

which is a generalisation of the Buchdahl solution (102) for b ^ 0. 

Deng in 1989 invented a powerful method (Deng's ladder) [3],[l5] for 
generating an infinite chain of more and more complicated solutions, varying 
with time, from a simple seed A\. One finds next the general form of L 1; 
takes it as a seed, finds the general A 2 and so on. He delivered the most 
general conformally flat solution and some others. Banerjee et al [61] gave 
two dynamical solutions 

T (t) z 1/2 - a z 

Al = — -H TTn , A H = — . . — . u , Z = l + 77(t)u. (Ill) 

T(t)^/2 + a' T(t) (z + a/T(t))' /W V 7 

In 1990 Knutsen [16] corrected and extended the Goldman potential 
method by studying the characteristics, the physical plausibility and the 
dynamical stability of the static models in terms of the Goldman- Knutsen 
generating function. He discussed the third Goldman solution and proposed 
a simpler one with g = au + b. 

Another static solution was given by Burlankov [62] 



L2= /„ + t -W' 

\u + b + V3/2j 
Recently, Pant et al [63] studied in detail the static metric 



(112) 



cos \/2b (d — u) „ , , , , , 

A= —, , , , B = cosec 2 (a + bu . 113 

sm (a + bu) v 7 v 7 



Finally, Msomi, Govinder and Maharaj [10j . with the help of Lie symmetry 
analysis, found five transformations, leading to new solutions from old ones 
and introduced their generating function W. 
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9 Conclusions 



We have given a global view upon the study of collapsing shear-free per- 
fect fluid spheres with heat flow. The application of the LG formalism has 
been advocated throughout the present paper. It provides a very compact 
formulation of the isotropy condition (12), namely Eq (29), and a very sim- 
ple expression for \&2 - Eq (56). The formulas for the other characteristics, 
Eqs (52-59), are also straight and tractable. Eq (56) clearly shows why the 
condition for conformal flatness is so similar to the isotropy condition. 

The LG formalism also presents the simplest possible version of the junc- 
tion condition. This has been demonstrated explicitly for conformally flat 
and for geodesic solutions. It gives the right functions to disentangle this 
condition into well known differential equations like the Abel equation, the 
Riccati equation, the Bernoulli equation or the linear one. This formalism 
yields an alternative derivation of the general solution when the metric com- 
ponents are functionally dependent. 

We have also discussed an unified study of separable solutions by incor- 
porating the simple linear in time ansatz into the general formula for the 
solution of the junction condition (76). 

One of the main objectives of the paper is to bring together the results 
of the static and dynamical group of authors, not only in the chronology of 
particular solutions, but in the discovery of generating functions. The recent 
proposition of the generating potential W [10] has prompted the search for 
similar functions, mainly in the work of the static group. A bunch of five 
generating potentials has been found, any of which provides the complete 
solution of Eq (29). Their common feature is the presence of the basic form 
AB. Its use reduces the LG equation to a first order or an algebraic one, 
depending on which potential of the pair is chosen as a seed. It seems that 
the Goldman-Knutsen generating function satisfies the simplest equation and 
a future task may be to continue the studies of Knutsen upon characteristics, 
pertinent to the dynamical models. Obviously, new generating potentials 
may be proposed by taking other functions of AB or its w-derivative. In view 
of this we hope that the enumeration of the existing generating functions, 
undertaken here, will prevent wasting of time in rediscoveries in the future. 

Finally, putting in order the four-dimensional case will help to lessen 
the efforts in investigating shear-free radiating collapse in higher dimensions 

mm- 
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